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We propose a generalized version of the Dielectric Breakdown Model (DBM) for generic breakdown 
processes. It interpolates between the standard DBM and its analog with quenched disorder, as a 
temperature like parameter is varied. The physics of other well known fractal growth phenomena 
as Invasion Percolation and the Eden model are also recovered for some particular parameter val- 
ues. The competition between different growing mechanisms leads to new non-trivial effects and 
allows us to better describe real growth phenomena. Detailed numerical and theoretical analysis are 
performed to study the interplay between the elementary mechanisms. In particular, we observe a 
continuously changing fractal dimension as temperature is varied, and report an evidence of a novel 
phase transition at zero temperature in absence of an external driving field; the temperature acts 
as a relevant parameter for the "self-organized" invasion percolation fixed point. This permits us to 
obtain new insight into the connections between self-organization and standard phase transitions. 
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Fractal growth phenomena, such as viscous fingering, 
electric discharges in dielectrics, fracture propagation 
and fluid flow in porous media, have attracted much at- 
tention in recent years. In the study of these phenom- 
ena many models - the Diffusion Limited Aggregation 
(DLA) |j, Dielectric Breakdown Model (DBM) @, and 
Invasion Percolation (IP) || to quote only but a few - 
have been introduced as a first step towards the under- 
standing of the dynamic emergence of fractal structures 
in nature. These models have been quite successful in re- 
producing the essence of the above complex phenomena 
by capturing some key ingredients. Their simplicity has 
also allowed for analytical studies of the scale invariant 
properties of corresponding structures j|,^J. 

To proceed further in this field there are two possi- 
ble directions: One is to identify possible interrelations 
between these models and understand their eventual uni- 
versality or fundamental differences. For example, these 
irreversible processes have a much lower degree of uni- 
versality than their counterparts in equilibrium statisti- 
cal mechanics, but still the basic "dynamical screening" 
mechanisms leading to fractal growth (as opposed to the 
growth of compact clusters) can be categorized in some 
generic classes, among them that arising from the exter- 
nal physical field (e.g. a Laplacian field in the DBM), or 
extremal dynamics and memory effects [^],[| . The second 
is to make these models more realistic or closer to real 
phenomena. For example, the interplay between these 
different mechanisms has been hardly studied so far. For 
instance, it has been recently found that the combination 
of Laplacian screening and extremal dynamics, results in 
a much stronger screening effect than that of each sin- 
gle mechanism, leading to fractal dimensions lower than 
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FIG. 1. Schematic diagram showing the limiting cases of 
our model as a function of temperature and the parameter n. 
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those characterizing the two different elementary mecha- 
nisms [0. 

Here we present a step forward in both of the previous 
directions. On one hand, we introduce a generalization of 
the available models in a broader class which permits us 
the identification of relevant parameters in more realistic 
growth phenomena. On the other hand this broader class 
is theoretically important as it is able to cast in a coher- 
ent framework apparently unrelated physical situations. 
This gives a new perspective that permits us to clarify 
the roles of the single physical mechanisms in competi- 
tion, namely external driving fields, quenched disorder, 
and thermal noise. The relative strengths of these effects 
are modulated in our model by three parameters, 77, a 
and T respectively. For some particular values of these 
parameters different well known growth models are re- 
covered (see figure Q). For T = we have QDBM, in 
the infinite T limit the model reduces to the DBM; at 
X] = and zero temperature we have IP, while for generic 
temperatures we have compact (EDEN ||) growth with 
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a temperature dependent correlation length. 

Let us now present the model, that we name as 
Generalized Dielectric Breakdown Model (GDBM). The 
medium in which a breakdown process propagates is dis- 
cretized as a regular (say hypercubic) lattice. The rel- 
evant degrees of freedom are placed on the bonds con- 
necting lattice sites. Each of these bonds can either be 
broken or not. Unbroken bonds are characterized by a 
given variable, Si, that can take different values (for in- 
stance, this can represent in a sketchy way the elonga- 
tion of a spring at i). The simplest situation one can 
think of is that of s, being a spin- like variable Si = ±1. 
We consider a cylindrical geometry with lateral periodic 
boundary conditions. This geometry, in fact, avoids all 
complications due to the discretized nature of the lattice, 
(like the anisotropy of growth along the diagonals of the 
lattice, relevant in geometries with radial symmetry H). 
We are interested in studying the process by which the 
bonds in the interface break down successively, and in 
how the cluster of broken bonds evolves. 

The dynamics proceeds as follows. First of all we as- 
sume that the breakdown process is quasi-static, i.e. only 
one bond is broken at each time step, and it is in the in- 
terface of previously unbroken bonds, this is, the cluster 
of broken bonds is connected. As initial condition we take 
as broken all the bonds in the lower row of the cylinder. 
Secondly, there are two distinct processes going on dur- 
ing the breakdown, the time scales of which are widely 
separated. The first, slow process, is the dynamics of 
breakdown itself. The second, fast process, is the relax- 
ation of the field configuration and of variables s$ after a 
breakdown event. In between two consecutives breakings 
the bond variables relax to their associated equilibrium 
state, whose associated energy or Hamiltonian will be de- 
fined later. The separation of time-scales is a common 
ingredient in many models for fractal growth and self- 
organized criticality || , and is physically a quite reason- 
able assumption. 

We now define the fast dynamics in between two con- 
secutive breakdowns. Every bond i is subject to a stress 
Laplacian field Ei. As in the DBM this (electric) field 
is given by the solution of the Laplace equation with the 
appropriate boundary conditions on the growing fracture 
[^],[| , namely the potential is 1 at the upper border of the 
cylinder (upper electrode) and on the lower one and on 
the broken bonds. This field acts over the bonds modu- 
lated by a parameter r\ E? , for r\ — any dependency 
on the external field is canceled. 

The second physical ingredient is quenched disorder. 
At each bond we define a random resistance Xi, dis- 
tributed as tt(x) — ai a_1 , x £ [0, 1] and a € [0, oo). This 
gives an idea of the bond tolerance to applied stresses; 
the larger Xi the larger the resistance. Other probability 
distributions can be used without affecting the properties 
of the model. We now define an effective local field ac- 
counting for the electric field and the disorder hi = / Xi 



which acts on the local degree of freedom generating 
an interaction energy given by 



H; — hi 



(1) 



There are looser bonds, with Xi large, which can accom- 
modate a larger stress E% induced by the Laplacian field, 
and tighter, or more fragile, bonds, with Xi small. 

In general one could also include nearest neighbor in- 
teractions, so that Hi — hiSi + J2j JsiSj to account for 
material rigidity. However, for the situations we are in- 
terested in, the local stress fields are much larger than 
the coupling J and J can be neglected. 

The equilibrium statistical mechanics of each bond 
variable among two successive breakdown processes is 
easily determined by introducing the temperature T. The 
partition function factorizes on the bonds: Z({hi},T) — 
U t z(h t ,T) and 
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-h iSi /T = 2 cosh(^/T). 



We can now calculate, as a function of hi and T, all 
averages. In particular 
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(si) = -Tq^: \nz(hi,T) 



taah(hi/T) 



(Ssj) = I-tanh 2 (^/T) 

are of interest. It is clear that (sj) is a measure of the 
stress exerted by the external field on the bond i and 
(6 s"i} measures the strength of thermal fluctuations. 

We consider the breaking probability Pi of a bond as 
an increasing function of the stress, and as a decreasing 
function of the amplitude of its thermal fluctuations; it 
is indeed when the local degree of freedom is thermally 
frozen in a stressed state, i.e. when it is very rigid, that 
breakdown is more likely to occur. On the other hand 
when the variance is large, the spin can better absorb 
stress and is more flexible. 

Guided by these considerations, the simplest dimen- 
sionless expression for the breaking probability Pj is 



Pi <x 



(fi) 



= svah(hi/T). 



(2) 



At each time step all the Pi are calculated, and with the 
corresponding probabilities one bond is broken. After 
that, the boundary conditions are automatically modi- 
fied and the field E updated at each point. Using the 
hypothesis of time scale separation, z(hi,T) and Pi are 
automatically updated, a new bond breaks down in the 
interface, and the process is iterated. 

The behavior of the model in some limiting cases can 
be analytically sorted out. 
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TABLE I. Fractal dimension of aggregates of size 
L = 64, 128 for different values of T (a = 1, n = 1). 



(i) For T ^S> 1 one can expand Eq. (||) around hi = 0, 
and find Pi oc E% /x,. In this limit the disorder does not 
play a relevant role and one obtains back the DBM (this 
same fact was already observed in ]Tc|]). 

(ii) For T < 1, on the other hand, the argument of the 
sinh is very large and essentially only the bond with the 
largest hi has a finite probability to break as T — > 0. In 
this way an extremal dynamics is generated, namely the 
quenched version of the DBM 

(iii) By setting rj — 0, the dependence on the physical 
field is canceled and we have a purely geometrical pro- 
cess. In fact, in that limit we obtain the IP extremal 
dynamics. For T > a characteristic length is gener- 
ated in the dynamics as we will show later by performing 
numerical simulations, and one observes asymptotically 
compact clusters, i.e. EDEN [0 growth. 

It is worth to stress that the previous limiting cases 
are quite general. For example one could also consider 
continuous degrees of freedom <jj taking values on a range 
[—Xi,Xi\. The specific equations change but the limiting 
form of P, given by Eq. (§) for T > 1 and T <C 1, 
remains unchanged. 

In order to study the properties of our model in more 
detail in generic cases, we have performed numerical sim- 
ulations. Each bond of the lattice is assigned a growth 
probability as in Eq. || and the dynamics starts from the 
lowest electrode, which represents the seed of the growth 
process. As soon as the fracture pattern reaches the up- 
per limit of the lattice, the dynamics is stopped, and a 
new realization starts. We performed a set of about 50 re- 
alizations of the dynamics of size LxAL with L = 64, 128, 
for a wide range of values of the temperature T, and dif- 
ferent values of 77 and a. For each set of realizations we 
compute by using a box-counting method |Q the frac- 
tal dimension, Df, of the resulting clusters. Data are 
collected only in the central L x 2L part of the lattice, 
which is sufficiently far from the lower transient regime 
and the upper interfacial region ||. In Fig. |^a-b we show 
a plot of Df(T) vs. T for L = 64 and L= 128 respectively 
and different values of a and 77, while in Tab. | we give 
some numerical values of Df(T) (for 77=1, a — 1). caca 

It can be verified that, as predicted, in the limits T — > 
and T — > 00 the model approaches the QDBM and 
DBM known values of Df, 1.15 and 1.66 respectively 
fifoU^j^]. The transition from the QDBM-like behavior 
to the DBM-like one is smooth and continuous (Fig. 
with continuously varying Df(T). We have checked that 



2.20 
2.10 
2.00 

u 

IX. 
H 1.70 
g- 1.60 
w 1.50 
1.40 
1.30 
1.20 
1.10 



2.20 
2.10 
2.0U 



1.50 
1.40 
1.30 
1.20 
1.10 





• L=64, T|=l, a=l 






-0 L=64, 71=1, a=0.5 




* 


-» L— 64 T|— S a—1 

. 3" " " " 


T T 


£ 


J I 

m ,-* 


1 : 

(a) 




10"' 10° 10 1 


10 2 10 




T 






• L=128,Ti=l ) a=l 






□ L=128,r|=l, a=0.5 






-a L=128,r|=0.5,a=l 






A 




- s 


. ..-.-.! " 


(b) : 



FIG. 2. a) fractal dimension £>/(T) of GDBM vs. the 
parameter T (dimensionless) for system size L = 64 (a) and 
L — 128 (b) with different values of 77, a. Note the irrelevance 
of a (disorder) for high temperatures. 
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FIG. 3. Results of the box-counting analysis for a very 
small temperature T = 0.02 and n — 0. The values 
D f = 1.97(2) for L = 64 and D f = 1.98(2) for L = 128 
we get correspond to asymptotically compact clusters and, 
are very far from the IP fractal dimension Df ~ 1.89. 



the critical properties of the model become less and less 
sensible to the value of a as the temperature T is raised 
(see Fig. 2), in agreement with our previous arguments. 
In our picture of the breakdown phenomena, this means 
that for very high temperatures, the thermal disorder is 
more relevant that the quenched one. 

When the Laplacian field is removed (77 = 0) the prop- 
erties of the model change dramatically. We observe a 
transition from an IP-like behavior to an EDEN like be- 
havior. For a value of T = 0.02 (the smallest value that 
we are able to implement in simulations), we get already 
compact growth, with Df — 1.98 ±0.02 for size L = 128; 
observe that the slight deviation from Df — 2 being a 
finite size effect) (see Fig. ||) This result is not affected 
by changes of the disorder strength. In fact, it can be 
easily shown that the extremal dynamics we recover in 
the T = limit, is independent on the value of a if there 
is no Laplacian field 

To analyze in more detail the crossover from IP to 
EDEN model we studied the correlation properties of 
the clusters in the following way. A fractal structure 
is characterized by the presence of voids of all sizes with 
no characteristic length; that is precisely the situation 
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for T = in our model. However, for T > voids are 
still present in the clusters but they have a characteris- 
tic, temperature dependent, size s c (T). Since voids are 
compact objects, the square root of s c (T) represents the 
correlation length £(T) of our model. We have studied 
the void distribution of our model for different small T 
and system size L = 128. By collapse plot techniques we 
have obtained the behavior of s c (T) = £ 2 (T) shown in 
Fig. |^, and we have found a power law divergence of the 
correlation length as we approach T = 0, with an asso- 
ciated anomalous exponent v = 1.07(6). This power law 
divergence of the correlation length can be interpreted as 
the hallmark of a second order phase transition along the 
temperature axis, with T c = 0. This is in line with re- 
cent results found for other self-organized models || . To 
test this hypothesis we studied the avalanche distribution 
D(s) for different T's (for a definition of avalanches see 
for example [Q]). This distribution is fitted by the scal- 
ing function: D(s) = s~ T f(s/ sq) where s is the avalanche 
size, r is the avalanche exponent, f{x) goes to a constant 
for x — > and falls exponentially for x — > oo, and sq is 
the cutoff size of the avalanches away from the critical 
point. By assuming s ~ T~~ ~ £ D t , and (s) ~ T~ 7 , 
where (s) = sD{s), we obtain the usual scaling rela- 



vD 



/. 7 



tions among avalanche exponents — 
Our simulations give (Fig.|): l/a = 1.69(5), 7 = 0.72(5), 
t = 1.60(2), Df = 1.87(3). These values together with 
the above estimation of v are in good agreement with the 
scaling relations among exponents, and support the in- 
terpretation of the transition at T = from IP to EDEN 
model as a second order dynamical phase transition (a 
more detailed report on numerical as well as analytical 
results will be given in a forthcoming paper. 

In summary, we have presented a general model which 
include DBM, QDBM, IP and EDEN models as limit- 
ing cases, and permits to interpolate among them by 
describing the the interplay between quenched disorder, 
thermal disorder and external stress fields. This repre- 
sents a step towards an unification of the models actually 
used to describe fractal growth phenomena, under a com- 
mon picture. We report a continuous change in the frac- 
tal dimension as temperature is risen interpolating from 
QDBM to DBM, and a novel phase transition from the 
extremal (self-organized) IP fixed point to EDEN type 
of growth. This can trace a connection between fractal 
growth dynamics and ordinary critical phenomena. 
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FIG. 4. (a): log 10 -log 10 plot of s c (T) = £ 2 (T) versus T for 
the GDBM model without field. We find £(T) ~ T~ v with 
v = 1.07(6). We interpret this as an evidence of a second 
order phase transition at T c = 0. (b): Scaling of so and (s) 
vs. T. 
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